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Abstract 

We investigate general thermodynamic stability conditions for the superfluid. This analysis is performed 
in an extended space of thermodynamic variables containing (along with the usual thermodynamic coordi- 
nates such as pressure and temperature) superfluid velocity and momentum density. The stability conditions 
lead to thermodynamic inequalities which replace the Landau superfluidity criterion at finite temperatures. 

1 Introduction 

Usually in experiments the vortices destroy superfluidity at velocities far below the Landau critical velocity. 
This is why the superfluid hydrodynamics equations can be expanded in powers of low velocities and one safely 
uses the first nontrivial terms of this expansion. 

Nevertheless, there is a number of experiments (see £Q) where the superfluid flow is investigated in small 
orifices. It has been shown that in these circumstances the maximum velocity is a decreasing function of 
the orifice width and may reach the order of the Landau critical velocity if the aperture is small enough. 
This means that all thermodynamic quantities of the superfluid become nontrivial functions of the not small 
superfluid velocity (i.e., it depends not only on the usual thermodynamic coordinates such as pressure and 
temperature). The only assumption one can make (and we do it) is that the fluid at rest is isotropic. This quite 
general statement of the problem is used in the paper; we find the complete set of thermodynamic inequalities 
in this light, i.e., the conditions imposed on thermodynamic functions for the superfluid to remain stable. 

Finally we employ the Landau phonon-roton model to calculate the highest velocity compatible with obtained 
thermodynamic inequalities and show that it can be interpreted as a critical velocity. This thermodynamic 
scenario supposedly explains the superfluidity break-up in small orifices. 

2 Stability 

When deriving general superfluid hydrodynamic equations it is usually supposed [2] that each infinitesimal 
volume of the liquid is (locally) in equilibrium and this equilibrium is stable. For the state of the liquid to be 
stable, it should provide an entropy maximum (at least local) for an isolated system. Instead of investigating 
the condition of the entropy maximality, it is convenient [2] to use another, equivalent to the first one, condition, 
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that is the condition of the energy minimality under constant entropy and additive integrals of motion. Thus, 
to examine if the state is stable or not, one must investigate the second variation of the energy. Such analysis 
will provide sufficient conditions for the energy minimality. 

Total energy of the superfluid E to t is an integral of the energy density E over the entire volume 

Etot = J Edv. (1) 

The energy density can be obtained via a Galilean transformation 

2 

E=^ + v s j + E . (2) 

Here v s is the superfluid velocity, p is the mass density and subscript denotes quantities measured in the 
frame of reference of the superfluid component (that is the frame where the superfluid velocity is zero). Namely, 
Eq and jo are the energy density and the momentum density (or, equally, the mass flux) with respect to the 
superfluid component. The former is a function of p, jo, and the entropy density S. Its differential can be 
written as 

dE = TdS + fidp + wdjo, (3) 

where Lagrange multipliers T, /x, and w are the temperature, the chemical potential, and the so-called relative 
velocity of normal and superfluid components. The liquid is isotropic and, consequently, the velocity w and the 
momentum density jo are parallel to each other, as expressed by 

w 

jo = j (T,p,w) — . 

w 

This leads to a useful identity for the partial derivatives of jo with respect to w: 



9 jo \ w k w l f dj \ ( 8 kl w k w l 
dw') Tp w 2 \dw) Tp 



Jo- (4) 



Further transforming J2J, we can rewrite it with the help of @ in the form 



dE = TdS + (a + 5s- - v s v n \ dp + (j - pv n ) dv s + v n dj, (5) 

where we denoted the total momentum density j = pv s + jo and the normal velocity v n = v s + w. 

As usual, stability implies that each "allowed" fluctuation increases the total energy of the system Etot- 
Allowed are the fluctuations leaving conserved quantities unchanged. This means that the minimality of E^ot 
must be investigated under fixed entropy and all additive integrals of motion: mass, momentum, and superfluid 
velocity. While the conservation of mass and momentum is well-known, conservation of the superfluid velocity 
worths a special comment. Really, since the superfluid flow is irrotational, the velocity v s is a gradient of a 
scalar: v s = V0. The same is true for the time derivative v s = V</>. This formula expresses the conservation of 
all three components of the vector 

V s = J v s dr. (6) 

Consider a macroscopic fluctuation of all the variables SS, dp, <5v s , and <5j. They are conserved and this 
ensures that the first variation of the total energy for a uniform system is identically zero 



SE io t= —) Sp+ —) 6v s + — *j dr = 0. (7) 
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The minimality criterion must be obtained as the condition of the positive definiteness of the second differential 
quadratic form. The matrix of this quadratic form is a Jacobian matrix 8x8: 



Q = 



d(T, v n , p + v^/2 -v s v n , j-pv n ) 



d(S, j, 



) 



Common rule states that it is positive definite if all principal minors Mi,M2, 
positive. We recursively test these minors: 

• The first positivity condition 

d(T, j, p, v s ) d(T, j, p, v.) d{T, v n , p, v s ) 



(8) 



. Ms in the top-left corner are 



Mi 



d(S, j, p, v s 



d(T, v n , p, v s ) (9(5", j, p, v s ) 



<9w 



as 



9jo 

<9w 



djo 
dT 



> 



corresponds to the usual requirement of the heat capacity positivity. It is shown below that (djo/dw) T > 
0, hence the last inequality eventually becomes 



dS_ 

dT 



p,tv 



djo 
dw 



djo 
dT 



> 0. 



p.w 



Positivity of the next group of minors is easily verified with the following transformation 



Q' = 



d(T, v n , p, v s 
d(S, j, p, v s 



d(T, j 1 p 1 v g 
d(S, j, p, v s 



d(T, v n , p, v s ) 
5(T, j, p,v s ) 



(9) 



(10) 



Whether the minors Mi, M3, M4 are positive is determined by the second multiplier in (|10|l . Required 
condition is therefore equivalent to the positive definiteness of the matrix 



di_ 



-1 




-1 















(dj /dw) Tp 
jo/w 

j /w 



Here we used ||1J| and chosen the direction of the w vector as the first coordinate. This adds to our 
collection two more inequalities 

jow > 0, 



djo 
dw 



> 0. 



(11) 
(12) 



The same transformation applied to the biggest minors gives: 



Q = Q' 



d(T, v n , (i + vl/2 - v s v n , j - pv n ) 



d(T, v n , 



) 



Q'Q' 1 
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Again, the minors Ms, Mq, M7, M§ correspond to nontrivial principal minors of Q" . We use the thermo- 
dynamic identity to relate the chemical potential p and the conventional pressure p 



dp = - - dT + v s dw. 
P P 

This gives 

a (p + z; s 2 /2 - v sVn ) \ = /dp\ = 1 /dp\ 

d P /T,v„,v a V^P/t : w P\ 9 p)t, w 

The following is an explicit representation of Q" sub-matrix corresponding to a four-dimensional space 
p, , v% , ; as before we let the a;-axis run along w direction. Using we obtain 



(dp/dp) Tw jp (dj /dp) T w - w 

(djo/dp) Tw -w p-(dj /dw) Tp 

p-jo/w 

p-jo/w 



Appropriate inequalities are: 

dp 

dp 



> 0, (13) 

T,w 



which is literally a generalized (to a non-zero inter-component velocity w) positive compressibility require- 
ment, 

jo < wp, (14) 

and 

dp\ / / dj \ \ f dj 



Inequalities ©, l|ll() . I|12|) . I|13|l . (|14|l . and l|15|l are sufficient conditions for the thermodynamic stability. 



3 Discussion 

In a "stopped- normal-component" arrangement, the mass flux f with respect to the normal component may 
become more convenient than jo — the mass flux relative to the superfluid one. The obvious relation between 
them / = pw — jo leads to the following reformulation of the inequalities: 

fw < 0, / < wp, (16) 



'dp 



< P 



(17) 
(18) 
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9 pJt,» \ dw Jp,T >P \ d P/w,T (19) 

As a simple application of the derived inequalities, consider them at w = 0. From (|16fl . (|17|l . I|18|l . and (|19fl 

we get 

!SL>* (IL- 

Using conventional notation, last inequality reads in the limit w — > 

p s > 0, p n > 0. (22) 



4 Phonon-Roton model 

Here we provide a usage example of the stability criteria for real superfluid 4 He. To calculate derivatives involved 
in the inequalities one take refuge in the microscopic approach. Simple and clear Landau phonon-roton model 
works pretty well in wide temperature and velocity ranges. We use this model to calculate the contribution of 
these quasiparticles to the "modified" free energy in the frame of reference of the superfluid component: 

F = Fo - wj . (23) 

Differential of this potential is given by 

dT = -SdT-jodw. (24) 
The modified free energy is obtained from the excitation spectrum with a conventional formula 

fi-'HM^))^ <25) 

We denoted the excitation energy e(p), which is given for two branches by the expressions 

e ph (p)=cp, £r (p)=A+ (p ~^° )2 . (26) 

Here and below, subscripts distinguish the quantities related to phonons and rotons, c is the sound velocity, A 
is the roton energy gap, m is the effective mass, and po is the momentum at the roton minimum 1 . A small 
dimensionless parameter mA/pQ ~ 0.03 1 ensures, e.g., that the Landau critical velocity is determined by 
v L = A/pq. 

!Data taken from 00: p = 0.145 g/cm 3 , A = 8.7K, m = 0.16m Hc , p = 3.673 10 s g~ 1 / 3 p 1 / 3 h, c = 23800 cm/s, dA/dp = 
-0.4710~ 14 cm 5 s- 2 , dm/dp = -0.45 10~ 23 cm 3 , dc/dp = 467 10 3 cm 4 s _1 g- 1 . 
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When integrated, these dispersion laws give the following contributions to the free energy: 

-2 



T 4 7T 2 f V? 

C 2 



— — 



90^ 3 c 3 

T 5/2 m l/2 

Po . wpo ( A 
— ; ; — r — smh — — cxp — — 

2 i/2 7r 3/2 ?i 3 w T y \ T 



(27) 
(28) 



One can obtain all 2 thermodynamic variables by differentiating this potential. Namely 

S = - 




Jo = - 

Inequality (|I5I) is the first to become invalid. Appropriate validity region is plotted in Fig. ^ The liquid is 
unstable above the curve. 



S3 




Figure 1: Critical velocity w c versus temperature T at normal pressure. Dashed line corresponds to the equation 
T = A — pow. Note that the condition T < A — pow holds true over entire stability domain. The "stability" 
critical velocity w c coincides with the Landau critical velocity ul at zero temperature and vanishes completely 
(|22Jl at the critical temperature T c (the A-point). In the phonon-roton model the critical temperature is T c w 
2.8 K. 

At zero temperature the critical velocity becomes the Landau critical velocity v^. It should also be noted 
that for systems where all quasiparticles can be described hydrodynamically (in other words, systems lacking 



roton branch) inequality l|15|) at zero temperature includes (dp/ dp)^ 



> i.e., w < c. 



J We neglect the quasiparticle contribution to the pressure derivative because it is just a small correction to the speed of sound. 
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5 Conclusion 



Experimentally, the superfluidity break-up in small orifices is believed to have the following nature (see £Q). 
Until the aperture size is too small the critical velocity does not depend on the temperature and increases as 
the size decreases. This is the very behaviour that is specific to the vortex-related critical velocity. 

When the orifice width is narrow enough the vortex-related critical velocity becomes so high, that the 
break-up scenario and its features change. The critical velocity does not depend on the aperture any more but 
decreases when the temperature increases. This behaviour is commonly associated (see [Q) with the Iordanski- 
Langcr-Fisher mechanism (see |B]). Nevertheless, this association lacks numerical comparison because no reliable 
information about the actual orifice shape is available. 

On the other side experimentally observed behaviour of the critical velocity can be attributed to the suggested 
stability criterion. In other words we provide an alternative explanation of experimental results based on an 
assumption that in narrow orifices the thermodynamic limit of w c is reached. 

We should also note that our approach to the critical velocity as a stability limit is similar to that used by 
Kramer[7|. Actually the inequality he employed is not a thermodynamic one. Moreover, generally speaking it is 
wrong. But numerical results for the critical velocity he obtained using the phonon-roton model do not deviate 
much from those plotted in Fig^ 
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